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ABSTRACT 



SN 1961V, one of Zwicky 's defining Type V supernovae (SN), was a peculiar transient in 
NGC 1058 that has variously been categorized as either a true core collapse SN leaving a black 
hole (BH) or neutron star (NS) remnant, or an eruption of a luminous blue variable (LBV) star. 
The former case is suggested by its association with a decaying non-thermal radio source, while 
the latter is suggested by its peculiar transient light curve and its low initial expansion velocities. 
The crucial difference is that the star survives a transient eruption but not an SN. All stars 
identified as possible survivors are significantly fainter, L opt ~ 10 5 L Q , than the L opt ~ 3 x 10 6 L Q 
progenitor star at optical wavelengths. While this can be explained by dust absorption in a shell 
of material ejected during the transient, the survivor must then be present as a L/« ~ 3 x 10 6 L Q 
mid-infrared source. Using archival Spitzer observations of the region, we show that such a 
luminous mid-lR source is not present. The brightest source of dust emission is only Ljr ~ 
10 5 L© and does not correspond to the previously identified candidates for the surviving star. 
The dust cannot be made sufficiently distant and cold to avoid detection unless the ejection 
energy, mass and velocity scales are those of a SN or greater. We conclude that SN 1961V 
was a peculiar, but real, supernova. Its peculiarities are probably due to enhanced mass loss 
just prior to the SN, followed by the interactions of the SN blast wave with this ejecta. This 
adds to the evidence that there is a population of SN progenitors that have major mass loss 
episodes shortly before core collapse. The progenitor is a low metallicity, ~ 1/3 solar, high 
mass, Mzams <^ 80M Q , star, which means either that BH formation can be accompanied by an 
SN or that surprisingly high mass stars can form a NS. We also report on the mid-IR properties 
of the two other SN in NGC 1058, SN 1969L and SN 2007gr. 

Subject headings: supernovae: general, supernovae: individual: SN 1961V, SN 1969L, SN2007gr 



1. Introduction 



We know that stars both explode, as core-collapse SN, and erupt in luminous transients that eject mass 
but do not destroy the star. In some cases, both types of transients produce similar, Type Iln spect ra, where 



the "n" indicates that the emission lines are narrow 2000 km/s) compared to a normal supernova (Schlegel 
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199CI iFilippenkd 1 1 9971) . Type Iln SNe seem to be cases where the blast wave is interacting with a dense 
circumstellar med i um created either by a massive wind or mass ejected in a pre-SN eruption (see, e.g., 
Smith et all (12008b . iGal-Yam et all (120071) and references therein). The mechanism of the erupt ions from 
LBV stars is not well-understood (see, e.g., iHumphrevs & DavidsorJ[l994l ISmith & Owocklbood) . but they 
eject material at velocities lower than normal SN. Unfortunately, the luminosities of the faintest SN are not 
well separated from those of the brightest eruptions, making it difficult to safely classify t ransients at the 



bound ary. These brightest of stellar eruptions are frequently referred to as SN "impostors" (|Van Dyk et al. 



2002). Correct classifications are important for understanding the rates and mechanisms of both processes. 
In p articular, we note t he recent debates about the nature of SN 2008S and the 2008 transient in NGC 300 
(see lPrieto et al.1 (120101) and references therein). 



The most obvious difference between the two cases is that the star survives only in the eruption sce- 
nario. Thus, there have been att empts to identify the surviving star for a number of the impostors, with can- 
didates identified for SN 1954J dsmith e t alhoOlLlv an Dvk et alJbooll). SN 1961V (see below) SN 1997bs 



Jvan Dvk et all 19991 . Li et alJl2002h . andSN 2000ch Jwagner et all2004l 



Pastorello et al. 



20101) . It is prob- 



ably safe to say that none of these identifications besides SN 2000ch is certain. There is, however, a second 
test. Most of the candidate survivors are fainter than the progenitors, and this is expected because the 
surviving star lies inside a shell of ejected material that probably forms dust as it cools. For the spectac- 
ular Galactic example o f r? Ca rina, ~ 90% of the emission is absorbed and reradiated in the mid-IR (see 
Humphreys & Davidsonl (11994)). Thus, a good test for these identifications is to find the mid-IR emission 
from the survivor and check that it matches the absorption indicated by the difference between the luminosi- 
ties of the progenitor and the survivorQ While some SN may be late time IR sources, they should evolve 
more rapidly and are unlikely to show the balance between progenitor luminosity, optical absorption and 
mid-IR emission expected for a surviving star. While frequently noted, this test never seems to have been 
carried out. We do so here for SN 1961V. 

The progenitor of SN 1961V was (likely) th e brightest star in NGC 1058, with m pg ~ 18 in the decades 



before the transient (IBertolall 19641 . IZwick y 1964). lUtrobinl (119871) estimated magnitudes in December 1954 
of B = 18.2±0.1, V = 17.7 ± 0.3, a nd B-V = 0.6 ±0.3. G iven a distance of 9.3 Mpc, the Cepheid distance 
to fellow group me mber NGC 925 (ISilbermann et al.ll 19961 ). and Galactic extinction of E(B- V) = 0.06 mag 
( Schlegel et all 1998 ). this corresponds to Mg ~ -12, making the progenitor one of the brig h test stars in any 



1964), 



(1199%. based on the data obtained by 



galaxy . Detailed discussions of the lig ht cu rve are presented in iDqggett & Branchl (|1985l) . iGoodrich et al 
1989i) . Humphreys & Davidsonl J 19941) . and|H umphre vs et al 



Zwickv 



Bertolal Jl963b . iBertolal d 1964b . iBertolal Jl967l) . bertola & Arpl Jl97pl) . and|FesenJ Jl985l) . Sometime 



between 1955 and 1960 the star started to brighten, reaching a plateau of m pg ~ 14 by the summer of 1961 
before briefly peaking at m pg ~ 12.5 in December 1961. At this peak, it was brighter than the maximum 
of the Type IIP SN 1969L (e.g. Idatti et al.|[l97lh and comparable to the peak of the Type Ic SN 2007gr 
(IValenti et al.ll2008b . the two other SN in NGC 1058. It then dropped in brightness, going through a series 
of extended plateaus, including a 4 year period from 1963 to 1967 where it was only moderately fainter than 



There can be problems in this accounting from binary companions (see lKochanek and chance coincidences 
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the progenitor, with m pg 19. After 1968 it had faded below the point of visibility, m pg <; 22. Spectra of the 
event were also peculiar (IBranch & Greensteinlll97lh . with relatively narrow lines (FWHM ~ 2000 km/s) , 
strong helium emission lines a nd a con s tant co lor, resembling an F star, near maximum light Jb ertolall 1 965 ) . 
Based on these peculiarities, IZwickyl (119641) classified SN 1961V, along with 7] Carina, as a "Type V" 
supernova, although th e clear presence of hydrogen in the spectr a would lead to a "modern" classification 
of Type Iln or peculiar (|Branch & Cowadll985UFilippenkdll997l) . 



Goodrich et all (I1989T) proposed that the peculiar, extended light curve and low velocity spectra would 



be more easily explained if SN 1961V was actually an LBV eruption rather than an SN. They proposed that 
the progenitor was a hot (T* > 45000 K), luminous (L* ~ 10 6 4 L Q ) star that was undergoing an S Doradus 
outburst in the decades prior to the eruption. During such an outbur st, the star has the same bolome tric lu- 
minosity but a far lower photospheric temperature (T* ~ 8000 K, see Humphreys & Davidson! (| 19941V ). This 
scenario requires a surviving star, and several candidates have been i d entified from a seque nce of steadily 



impr oving Hubble Space Telescope i mages b y Filipp enko et alj <jl9 95), Van Dv k et all ( 20021) and 



(2004 1) based on the a ccurate optical (|Klemolalll986T) or radio (|Branch & Cowan 



1985 



Chu et al. 



Cowan et al. 



1988, 



Stockdale et al.ll200lh positions. All proposed candidates are significantly fainter than the progenitor, with 
V ~ 24 mag. 

The primary counterargument to the LBV eruption hypothesis is that SN 1961V also seems to be asso 



dated with a fading, n on-thermal radio source (|B ranch & Cowanll 19851. ICowan et al.lll988L IStockdale et al. 



200 ll. IChu et al.ll2004f ) that closely resembles the properties of other radio SNe and not the fainter, thermal 
emission of rj Carina. There is no non-thermal radio emission (or even a detection) fro m the SN impos - 
tor/LBV eruption SN 1954J ev en though i t is alm ost four times closer and of similar age (|Eck et al.l l2002). 
VLBI observations in 1999 by IChu et alj (120041) also resolved out the radio emission, setting a minimum 
radius for the radio emission of order 4 mas or about 0.17 pc. This would require an expansion velocity of 
v <; 4000 km/s that would be hard to explain with an eruption. 

While the Spitzer Space Telescope (SST) was not intended for studies of individual stars at 10 Mpc, 
it should have no difficulty ident ifying a source with the luminosity of the SN 1961V progenitor star in the 
outskirts of NGC 1058. Indeed, boodrich et aP dl989h note that in the infrared the source should be "the 
brighte st point thermal I R source in NGC 1058." Here we use archival SST IRAC qFazio et all 120041) and 
MIPS (IRieke et al.ll2004l) data to measure the infrared emission associated with SN 1961V. In $2]we discuss 
the available data, the astrometry relative to the HST data used to identify candidate surviving stars, and the 
resulting estimates and limits on the mid-IR luminosity. We model the photometry in $3] to find that there 
is insufficient infrared emission for the progenitor s tar to have survived and that SN 1961V must, therefore 
have been an SN. We note that lSmith et al. ( 2010 ) have simultaneously reached this conclusion based on 
the gross differences between SN 1961V and other supernova impostor candidates and its greater similarity 
to other core-collapse SN. In §@]we present the photometry for SN 1969L and SN 2007gr, the other two SN 
in NGC 1058. In §f5]we discuss the consequences of SN 1961V having been an SN. 
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Fig. 1. — Astrometric matches between the F606W HST image and the 3.6/xm reference image. The 
top panels show a 45" field of view showing three stars (labeled A, B and C and marked by I'.'O radius 
circles) that can be well-matched between the bands. The lower panels show a narrower lOf'O region 
around SN 1961V co rresponding to the box in the top panels, where we have labeled the so urces following 
Van Dyk et all (120021) . The region labeled X marks the location of stars #6, 7, 9 and 1 1 from lVan Dyk et al. 
(120021) and contains all the candidate surviving stars. Source #6 is brighter in the F814W image, and sources 
#3 and 11 are only detected in the F814W image. The small circles in the lower panels are 0"3 in radius 
while the large circle corresponds to the 2"4 radius photometry apertures shown in Fig. [2] 
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Fig. 2. — Mid-IR images and wavelength differenced images of the SN 1961V region. The top panels show 
the 3.6, 4.5 and 5.8/im images of the region, the middle panels show the 8.0, 24 and 70/mi images of the 
region, and the lower panels show the [4.5] -[3.6], [5.8] -[3.6] and [8.0] -[3.6] wavelength differenced 
images. This removes the flux from normal stars to leave only sources of dust and PAH emission. The 
panels are 15" 0, 30"0 and 60f in size for the IRAC, 24/im and 70/xm bands, respectively. The large black 
circles in the IRAC, 24 and 70/t/m panels have radii of 2'.'4, 3'.'4 and 16"0, respectively, and correspond to 
the aperture sizes used for photometry. The smaller 1 '.'2 radius circles mark the positions of star #8 and 
the region X encompassing the candidate surviving stars. We also analyzed the IRAC images of the region 
using a 3 "6 radius aperture and DAOPHOT. 
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Fig. 3. — Mid-IR SEDs of the SN 1961V region. The total emission is described by the large aperture 
MIPS fluxes and either the 3 "6 (open squares) or 2 "4 (filled squares) IRAC apertures. With DAOPHOT 
we attempt to separate the fluxes of star #8 (filled triangles) and the region X (open triangles) that contains 
all the proposed surviving stars . For comparison, we show with open pentagons the SED of rj Carina from 



Humphreys & Davidson! (119941) . which roughly has the properties we expect for SN 1961V. The 2'.'4 IRAC 



apertures combined with the 24/mi luminosity and the 70/mi upper limit will be our standard comparison 
SED. Note that region X has an SED dropping to longer wavelengths, indicating it is dominated by stellar 
emission, while star #8 has an IR excess. 
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2. Data and Infrared Luminosity Estimates 



NGC 1058 has been observed twice with IRAC and MIPS, as summarized in Table [T] The total ex- 
posure times are 15 x 30 s = 450 s for the IRAC bands, 10 s + 30 s = 40 s for the MIPS 24 fim band and 
3 x 3 s = 9 s for the MIPS 70/mi band. The SST sensitivity estimates for these exposure times are 0.36, 
0.62, 4.1, 4.7, 25 and 3400^ Jy, at 3.6, 4.5, 5.8, 8.0, 24 and 70/mi, respectively. If we convert these into 
3a limits on uL v in each band at the distance to NGC 1058, they correspond to 2400, 3300, 17000, 14000, 
25000 and 1.2 x 1O 6 L for the 3.6, 4.5, 5.8, 8.0, 24 and 70/mi bands, respectively. In practice , we would be 



confu sion limited in the IRAC bands were we trying to reach these detection limits, but, as iGoodrich et al. 



(1989) noted, we should have little difficulty finding the expected > 10 6 L Q mid-infrared source! 



We downloaded the Post-Basic Calibrated Data (PBCD) for these programs from the Spitzer archive. 
These IRAC images are two-times oversampled and have a pixel scale of 0"60, while the MIPS 24/mi and 
70/mi images have pixel scales of 2"45 and 4"0 respectively compared to native pixel sc ales of 2"55 and 
5 "2 (narrow field o f view). We aligned and combined the data for each band using the ISIS (lAlard & Lupton 



1998 



Alardll2000|) image subtraction package. We also used ISIS to difference image between the avail- 



able epochs, to search for any signs of variability, and to difference image between wavelengths. The latter 
technique takes advantage of the fact that all "normal" stars have the "same" mid-IR colors, so normal 
stars effectively "vanis h" to leave only t he red stars dominated by dust emission and emission by the in- 
terstellar medium (see IKhan et al.ll2010i) . This wavelength differencing procedure isolates the relatively 
rare, dusty stars with out the crowding from the normal stars. We also obtained the HST images used by 
Van Dvk et al.l (120021) so th a t we could astrometric ally match the Spitzer data with the progenitors discussed 
by iFilippenko et all dl995h . Ivan Dyk etaD (I2OO2I) and bhu et all d2004l) . We also examined the more re- 
cent images of the area from October 2007 (Van Dyk/1 1119) and August 2008 (Li/10877), but these do not 
significantly improve on the prior data. 

Fig. [I] shows a wide field and close-up view of the SN 1961V region in the HST WFPC2 F606W 
(Illingworth/5446) and 3.6^m reference images. Clearly, with Spitzer's resolution we will be unable to 
obtain photometry for all the individual stars identified in the HST image, particularly at the longer wave- 
lengths. We see counterparts in the 3.6^m image to star #8, the group of stars #5/6/7/9/11 (which we will 
refer to as region X) and star #3. Stars #3 and #11 are not visible in the F606W image, but are detected in 
the F814W image. The large black circle is 2"4 in radius and represents one of the apertures we used for 
photometry. 

Fig.|2]shows regions around SN 1961V for all 6 Spitzer bands. We used black circles to mark the 2"4 
photometric aperture we used to estimate the fluxes on IRAC images, as well as the 3 "5 and 16" apertures 
used to measure fluxes in the 24.0 and 70.0/im MIPS bands respectively. For clarity, the alternative 3"6 
aperture used for the IRAC bands is not shown. We also show the wavelength differenced images between 
3.6/um and the other three IRAC bands. We see that most of the sources in the 3.6/mi image are normal 
stars, since they fade away at longer wavelengths and do not appear in the wavelength differenced images. 
There is some dust related emission, much of which seems to be associated with source #8, for which we 
lack an optical color because it lay just off the field edge in the F814W and F450W images analyzed by 
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Van Dyk et all (120020 . and possibly with source #10. IChu et al.1 (120041) identify star #7 as the only point-like 



source of Ha emission. The complex of sources corresponding to stars #6, 7, 9 and 11 in IVan Dyk et al. 



(2002), which we have labeled region X in Fig.Q] appears to have no significant excess emission due to dust 
even though they correspond to all the claimed counterparts to SN 1961V (see below). 

We estimated the fluxes using two procedures. First, we simply used aperture photometry (the IRAF 
apphot package). We used signal aperture radii (background annuli) of 2"4 (I'.'A-l'.'l) and 3"6 (3"6-8"3) 
for the IRAC bands, 3"5 (6f'0-8f'0) at 24/mi and 16"0 (18;'0-39f'0) at 70/xm. The background was esti- 
mated using the mode of the background pixels after 2a outlier rejection, an approach which should work 
reasonably well in crowded regions. We also compensated for the presence of the edge of the 70/xm image. 
No source was identified at 70/mi, so we estimated a 3er upper limit on the flux. We used the standard 
Spitzer corrections for these apertures]! For the 2'.' A (3" 6) IRAC aperture these are 1.213, 1.234, 1.379 and 
1.584 (1.124, 1.127, 1.143 and 1.234) for the 3.6/im, 4.5/xm, 5.8/im and 8.0/um bands, respectively, with 
uncertainties of order 1-2%. For the 24 and 70/xm apertures, they are 2.80 and 2.07 and are accurate to 
about 5%. The resulting flux estimates are presented in Tables |2] and [3] 

While the large aperture photometry provides a conservative upper limit on the luminosity of any 
individual source, it is clear that the flux near SN 1961V can be divided over several sources in the IRAC 
images. To better account for the effects of overlapping PS Fs than is possible with aperture photometry, 
we also analyzed the region with DAOPHOT (|Stetsonll 19871 ). in particular dividing the IRAC flux between 
source #8 and the complex of sources in region X associated with the candidate surviving stars. For the 
still lower resolution MIPS images, no attempt was made to divide the flux over sub-components. The 
DAOPHOT results are also presented in Table 12 Fig . |3] com pares these estimates to each other as well as 
to the SED of t] Carina from lHumphreys & DavidsonI (119941) . The total flux, even in the large apertures, is 
far less than that of i] Carina, which roughly has the luminosity and SED we expect for SN 1961V. The sub- 
components are then significantly less luminous, although the DAOPHOT division into two sources does 
not capture all the flux in the aperture, and we again see that while star #8 has an IR excess, the region X 
containing all the proposed surviving stars seems not to. We also checked for variability in the IRAC and 
24/xm bands, finding none to limits of roughly 10%. 



We also repo rt photometry f or Stars A, B and C in Fig. [Tl IVan Dvk et all (|2002h star #3, SN 1969L 
(ICiattietal.ll 197 lb and SN 2007gr Jcrockett et allboOsllyaTenti et al.ll2008h . The location of SN 1969L was 
only covered by some of the images and SN 2007gr is only present in the later data, so we only analyzed 
the relevant images but followed the same procedures. We only obtained upper limits on any flux from 
SN 1969L, while SN 2007gr was a very bright source. We discuss the results for these SN in $4] 



2 http://ssc. spitzer.caltech.edu/irac/iracinstrumenthandbook/ and http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/ 
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3. Models 



We model the spectral e nergy distributions (SED) using DUSTY dlvezic&Elitzur!ll997l . Ilvezic et all 



19991 . Elitzur & Ivezid 120011) . We assumed a dusty shell with a density distribution oc 1/r 2 and an outer 
radius at twice the distance of the inner, R out = 2Ri n . This assumption has little affect on the results. The 
models are specified by the temperature of the illuminating black body, the stellar temperature T*, the optical 
depth of th e shell Ty at V band, and the dust temperature Tj at the inner edge of the shell. We tabulated the 
models for lDraine & Led (| 19841) graphitic and silicate dusts with the standard size distributions assumed by 
DUSTY for stellar temperatures of = 5000, 7500, 10000, 15000, 20000, 30000 and 40000 K, inner edge 
dust temperatures from 50 to 900 K in steps of 50 K, and V-band optical depths of t v = to 6 in steps of 0.1 
and ry = 6 to 30 in steps of 0.5. Our approach will be to normalize the models based on the pre-transient 
luminosity and then constrain the optical depth to match the flux of the candidate survivors, leaving as the 
remaining variable the dust temperature. 

Given the stellar luminosity, the dust temperature is determined by the shell radius, and the shell radius 
is closely related to the physics of the transient. Since the progenitor must have been essentially unobscured 
pre-transienH, we know that any obscuring material must have bee n ejected in 1961. If the v elocity is 
restricted by the FWHM of the optical lines, roughly 2000 km/s (e.g. Branch & Greensteinlll971l) . then the 
current radius of the material is 



R~ 1.4 x 10 



17 



cm 



(1) 



1000 km/s 

where we set the elapsed time to 43 years (1961 to 2004) and the velocity v e j to half of the FWHM. In our 
standard models we consider those with inner shell radii near this value, which is mildly conservative given 
that the inner edge dominates the dust emission. The outer edge is then at twice this distance and so has 
twice the expansion velocity. 

The only way to escape our eventual limits is to make the dust so cold that it cannot be detected given 
Spitzer's diminishing sensitivity at longer wavelengths. For the dust temperature to be low, the dust must 
be distant, and for the simple case of radiative equilibrium for dust radiating as a black body, the dust 
temperature is 

/ L \ 1/4 f L \ 1/4 / 10 17 
r= (_^_) = i42 / ' N /W 



3 x 10 6 L m 



cm 



R 



1/2 



K 



(2) 



(e.g. IWrightl (119801) ) corresponding to an SED peaking near A = 20fim that will be strongly constrained 
by the 24/xm data. DUSTY, with better dust emissivity models, usually predicts higher inner edge dust 
temperatures than this simple model but a similar peak emission wavelength. Lowering the dust temperature 
to raise the peak wavelength requires a larger dust radius, but moving the shell to a larger radius requires 



3 Otherwise, its already high luminosity would quickly exceed 10 7 Lq after extinction corrections! Adding additional unrelated 
foreground extinction only strengthens our conclusions because it will adjust the progenitor luminosity upwards without contribut- 
ing to the mid-IR flux. 
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rapid increases in both the ejected mass and energy. The V-band optical depth of the shell is 



T V ■ 



M e jK opt 

4irR ir ,R im1 



^opt 



500cm 2 /g7 V Ri< 



Rn 



10 17 cm 

Rin 



(3) 



where K opt ~ 



k 500 = 500 cm 2 /g is the optical opacity for a dust to gas ratio of roughly 1% (e.g. Isemenov et al. 



2003), M e j is the ejected mass in which the dust forms, and the shell has a density profile oc 1 /r 2 from 
Ri n <R< R ou t- Equivalently, the required mass is 



M £7 = 0.13r y 



^500 



^opt 



Rn 



Rn 



R„ 



10 17 cm 



M G . 



(4) 



Assuming a thin shell with /?,•„ ~ R out = R for simplicity, the energy of the ejecta, 



-M ej v 2 ej = 7 x 10 47 tv 



^500 



^opt 



R 



10 17 cm 



ergs, 



(5) 



increases very rapidly with increasing shell radius because both the velocity and mass must be larger for 



larger distances. The energy required reaches an SN-like magnitude of 10 M ergs for R ~ 6 x lO 17 ^ 1 ^ 4 cm, 
as does the velocity and mass. Making the dust distant enough to be cold forces the mass, velocity and energy 
budgets out of the LBV transient range. Phrasing the scaling in terms of the peak wavelength, M e j oc \ A peak 
and E e j oc \ 8 peak , further emphasizes the problem with this solution. Using a thick shell exacerbates these 
problems since it leads to a larger mass-weighted radius. 



In order to understand the mid-infrared limits we must first choose which HST star to call the survivor 



(see Fig. [T])- Filippenko et al 



995|) propose star #6, with B and V magnitudes of 24.82 ± 0.25 and 24.50 ± 



0. 16 mag. IVan Dyk et al.1 (|2002 ) propose s tar #1 1, which they estimate t o have an I magnitud e of 24.3 ±0.22, 



B-I > 1 and V-I > 1.1 mag. IChu et all (|2004h propose star #7 from I Van Dvk et all (120021). wi th B, V and 
I magnitudes of 24.04 ±0.14, 23.85 ±0.14 and 23.83 ±0.14, respectively. Ichu et all J2004h prefer this 
identification because (1) it appears to be spatially coincident with the radio source, (2) it appears to be 
the only Ha source, and (3) its Ha line has broad wings (at least ±550 km/s, limited by the noise in the 
spectrum). It does not, however, have the forbidden lines ([OI] A6300A, [OIII] A4959/5007A) expected 
from a remnant, sugge sting the Ha emis sion is st ellar. It is also too blu e to be well-modeled as an extincted 
(hot) star, as noted by IChu et all (12004) . Of the IVan Dyk et al.l ((2002) candidates closest to the preferred 
positions, #5 and #9 also have the wrong spectral slopes, while #6 and #11 are easily fit. Star #8, which has 
not been proposed as a candidate because it is too distant, is the only nearby source with significant dust 
emission (see Figs.[2]and[3]>. 

In practice, it matters little which star we use for the survivor. They are all faint compared to the 
progenitor star and so must be heavily extincted in the optical. Once most of the optical/UV flux must 
be absorbed, it matters little for the expected mid-IR luminosities whether it is 90% or 99%. Similarly, 
the spectral shape of the candidate matters little, since the optical depth is principally determined by the 
magnitude difference between the survivor and the progenitor rather than the color. Thus, for simplicity we 
will simply give the survivor the typical magnitude of the candidates, V = 24 mag, and ignore the colors. 
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Fig. 4. — Models normalized by the pre-outburst magnitudes of lUtrobinl (|1987h . The predicted SEDs for 
graphitic (silicate) dust shown by the heavy solid (dashed) curves lie far above the total mid-IR emission 
(filled squares) from the SN 1961V region let alone that of any sub-component (see Fig. [3]). In this case, 
the progenitor model (light solid curve) i s a = 7500 K black body normalized to match the pre-outburst 
magnitudes (open squares) from lUtrobinl (119871) . The progenitor luminosity is L* = 10 6 9 L Q and it would 
increase, leading to larger discrepancies, if we used a higher or lower stellar temperature. The V band optical 
depths are chosen to match the luminosity corresponding to our generic V = 24 mag extincted, surviving star 
(filled pentagon), and the inner shell radius is set to be close to 10 17 cm. 
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Fig. 5. — Models following the [Goodrich et alJ (|1989T ) scenario. The predicted SEDs for graphitic (silicate) 
dust shown by the heavy solid (dashed) curves still lie well above the total mid-IR emission (filled squares) 
from the SN 1961V region let alone that of any sub-component (see Fig. [3). In t his case, the progenitor 
model (light solid curve) is a = 30000 K black body normalized to match the [Goodrich et alJ (119891) 
normalizing magnitude of B = 22 mag, leading to a stellar luminosity of L* = 10 6 2 L that is somewhat 
low. The V band optical depths are again chosen to match the luminosity corresponding to our generic 
V = 24 mag extincted, surviving star (filled pentagon), and the inner shell radius is set to be close to 10 17 cm. 
For T* = 40000 K, L* would double and be closer to matching the luminosity in Fig. HJ which would also 
double the mid-IR discrepancy. 
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Fig. 6. — Stellar luminosities and shell radii (left scale) or ejecta velocities (right scale) that fit both the 
generic V = 24 mag extincted luminosity of the candidate survivors and roughly stay below the upper bound 
on the mid-IR emission (x 2 < 24 in Eqn.© for cold (T* = 7500 K, top) and hot (T* = 30000 K, bottom) stars 
and either graphitic (left) or silicate (right) dust. The vertical lines indicate the minimum (ry = 0) luminosity 
consistent with a V = 24 mag survivor. The open pentagon marks the solution with the properties typically 
associated with LBV transient hypothesis at L* = 10 6 2 Lq and v e j = 1000 km/s. 
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Fig. 7.— Limits on the mid-IR SED of SN 1969L (filled t riangl es) and the observed SED of SN 2007gr 
(filled squares). The filled R-band point from IValenti et alj (I2008I) corresponds to the epoch (+17 days) of 
the IRAC observations while the open point corresponds to the epoch (+3 1 days) of the MIPS observations. 
The heavy solid curve is a 5000 K black body with luminosity 10 8 4 L Q . 
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We have two possible choices for the intrinsic properties of t he star. First, we can simply normalize 
it using the pre-transient luminosities. Alternatively, we can follow iGoodrich et al. d 19891) and assume that 
the star was in an S Doradus phase just before the transient and has now returned to its hotter, bu t similar 
lumin osity quiescent state. If we set the stellar temperature to = 7500 K and nor malize it by the lUt robin 
(119871) magnitudes, we get a luminosity of L* ~ 1O 6 9 L . This is higher than the IGoodrich et al.1 (I1989T) 



proposal of 10 6 4 L Q , but matches their p roposed S Doradus outburst temperature and agrees tolerably well 
with the B-V color from lutrobin (1987). Both raising and lowering the assumed temperature increases the 
luminosity because for = 7500 K the peak of the SED lies near the normalizing photometric bands. In 
fact, with any significant change in the temperature, the luminosity becomes impossibly large (L* > 1O 7 L ). 

Fig. |4] shows the resulting models. Dust optical depths of ry = 7 and 11.5 for graphitic and silicate 
dusts lead to enough extinction to make the optical flux consistent with a V = 24 mag flux for the surviving 
star. Choosing inner edge dust temperatures of Td = 500 and 300 K corresponds to putting the inner edge 
for the dusty shell at R ~ 1.3 x 10 17 cm. The outer edge dust temperatures are roughly 140 K. Using a 
thin shell, with R out = 1.2/?,„ instead of 2R[ n leads to no significant changes. It is immediately apparent that 
the predicted mid-IR emission is grossly discrepant with the constraints even when we compare the model 
to the integrated emission from the region (defined here by the 2'! A IRAC fluxes and the MIPS aperture 
fluxes) without any division of the emission over the multiple sources within it (see Fig©. The cool stellar 
temperature exacerbates the problem because much of the near-IR emission is little affected by extinction. 



The alternate hypothesis, that the star has reverted to a quiescent hot state ( Goodrich et al.l ll989). 
changes things little. Here we assume that the surviving star has now left its S Doradus phase, and again has 
a high photospheric temperature with a quiescent magnitude that is 4 mag fainter than before the eruption, 
B ~ 22 mag. For black bodies with = 30000 K and 40000 K, this implies luminosities of L* = 1O 62 L 
an d L* = 10 6 5 L^ that are s ignificantly below that implied by the pre-eruption luminosities, as also noted 
by Humphreys et al.l (I1999T) . Fig. [5] presents the = 30000 K models. The visual optical depths are now 
much smaller (-ry = 2.5 and 4.5 for graphitic and silicate dusts) because most of the flux is in the UV where 
the dust opacities are higher. Choosing inner edge dust temperatures of Td = 400 and 300 K corresponds to 
putting the inner edge for the dusty shell at roughly R ~ 1.5 x 10 17 cm. The outer edge dust temperatures are 
roughly 100 K. Again, using a thin shell leads to no significant changes. The discrepancies are smaller here, 
partly because the stellar luminosity is a factor of 5 lower than in the models of Fig.|4j and partly because the 
star has little near-IR luminosity compared to the cooler model. Nonetheless, the predicted mid-IR fluxes 
are still much higher than allowed by the observations. Raising the stel l ar tem perature to 7* = 40000 K in 
order to better match the pre-transient luminosity or the IGoodrich et all (119891) models makes the problem 
worse by a factor of two. 

In the end, the two most important variables are the intrinsic luminosity and the radius of the dust 
shell, or equivalently the ejection velocity of the material. We can explore these models by normalizing our 
DUSTY models to fit the generic V = 24 mag of the obscured progenitor star and then keeping only the 
models consistent with the total luminosity of the region, L fo(A,), again defined by the 2"4 aperture IRAC 
luminosities and the MIPS aperture luminosities and limits. The shell radius, or equivalently the expansion 
velocity, is the primary secondary variable. We can quantify the consistency using the mid-IR luminosities 
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for the region as luminosity limits based on the metric 

X 2 = ^(L mod (A i -)/L ofo ,(A i )) 2 (6) 

i 

where A,- corresponds to the 6 Spitzer bands, L fo(A,) is the Spitzer luminosity limit and L mot i(Xi) is the 
luminosity predicted by the model. An SED passing exactly through these 6 values or limits would have 
X 2 = 6, while those that are too bright will have higher values and those that are fainter will have lower 
values. At least for the IRAC bands, we know that these large aperture luminosities are subdivided over 
multiple sources and so are upper limits even if any of the source identifications are correct. In Fig. [6] we 
show all cases with x 2 < 24, which corresponds to the typical model luminosity being twice the observed 
upper limit, for stellar temperatures of = 7500 and 30000 K and for graphitic and silicate dusts. As the 
secondary parameter we use the radius of the inner edge of the shell, or equivalently the expansion velocity 
of the shell. We only show cases with ry > 0.1. 

If the star is hot, T* = 30000 K, then the stellar luminosity has to be at least L* > 10 5 2 L Q in order to 
fit the optical luminosity of the surviving star. This is already high enough to make it very difficult to fall 
below the upper limits on the mid-IR luminosities in the presence of any dust. A few silicate models work by 
putting the shell so far out that the dust is cold enough to avoid all but the weak 70^m limit or by making the 
optical depth significantly less than unity. More solutions are possible when the star is = 7500 K, because 
the star is intrinsically less luminous for the same V band normalization. We see the expected trend allowing 
more luminous stars for more distant and colder shells. There are no solutions in the regime required by the 
LBV eruption hypothesis. 

We see no signs of variability in the mid-IR at the level of about 10% of the IRAC fluxes, although we 
would not expect to given the limited time baseline. While the existing data is not of high enough quality to 
make the test, we note that the optical variability should be significant given the parameters of our models. 
As the optical depth drops, the source should becom e steadily brighter, as is observed for rj Carina (e.g. 



Humphreys & DavidsodI 19941 Humphreys et al.lll999T) . If we normalize the optical depth to r\ and time t\, 



the optical depth scales as r = T\{t\/t) 2 (Eqn. [3]>, although the expected magnitude does not simply scale 
with r because it includes both scattering and absorption. If we take Ty = 4.5 from the silicate models for a 
hot star, then a V = 24.0 mag progenitor in 2004, then it should have been 25.7 mag in 1995, and should be 
23.4 mag in 2010, 23.1 mag in 2015, and 22.9 mag in 2020. In the cool star models the evolution is even 



more dramatic because of the higher optical depths. As IChu et all (120041) noted, there is no sign of optical 
variability in the published data. 



4. SN 1969L and SN 2007gr 

Since we were analyzing the Spitzer data already, we measured the fluxes associated with the other 
two SN in NGC 1058, SN 1969L and SN 2007gr, reporting their fluxes in Tables [2] and [3] and presenting 
their SEDs in Fig|7J SN 1969L lay outside the 70/im image, and the image was taken before SN 2007gr, 
so we have no information on their 70/xm fluxes. We detected no flux above background at the location 
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of SN 1969L JCiatti et al.lll97lh . at limits on vL v of 10 4 to 1O 5 L . Here we used the 3a upper bounds 
from the 3'.' 6 radius IRAC aperture and t he normal 24/mi aperture. The observations of Type Ic SN 2007 gr 
dCrockett et alJbood. Ivalenti et al.ll2008h were taken 17 (IRAC) and 31 (MIPS) days after the R-band peak. 
The mid-IR luminosities are comparable to what was expected for SN 1961V, with luminosities vL v of 
10 70 , 10 6 - 8 , 10 6 - 5 , 10 6 - 4 and 10 6 - 1 for the 3.6, 4.5, 5.8, 8.0 and 24/mi bands respectively. The SED, as 
shown in Fig. 13 is falling rapidly in the mid-IR, indicating that the emission is not dominated by a cool 
dust echo . We estimated R -band magnitudes at the two epochs of 13.6 and 14.3 mag based on the light 



curve m 



Valenti et all (120081) . The combined optical and mid-IR SED is well fit as a 5000 K black body 
with a luminosity of 1O 8 4 L . If we look at the wavelength differenced 4.5fim image from before the SN, 
there does appear to be excess emission, consistent with the presence of the stars having K band excesses 



in 



Crockett et all (120081) . but with the resolution of Spitzer and the presence of bright stars just North and 



South of the site, we cannot say more. 



5. Discussion 

The basic conundrum is simple. All possible surviving stars are far fainter than the progenitor in the 
optical. This requires significant visual optical depths so that most of the surviving star's luminosity is re- 
radiated in the mid-IR. However, this opacity must be supplied by the material ejected during the transient, 
and dusty shells ejected at the low velocities of the LBV hypothesis and irradiated by a surviving star lead to 
mid-IR luminosities in gross conflict with the observational limits. The discrepancy is not a subtle problem, 
but a disagreement of an order of magnitude or more. The limits could be evaded by pushing the shell so 
far outward in radius that the dust temperature is too low to have significant emission in the 3.6 to 24/^m 
range, but this solution requires shell masses, velocities and energies that are extreme even for a true SN. 
The simplest solution to these problems is that SN 1961V was in fact a supernova and there is no surviving 
star. 

To escape the conclusion that SN 1961V was a supernova requires that an assumption about the prop- 
erties of the star or its surrounding dust is greatly in error. We can enumerate three possibilities for changes 
in the stellar properties: (1) the system was (bolometrically) super-luminous for ~ 3 decades prior to the 
transient; (2) the surviving star has been (bolometrically) sub-luminous for the ~ 5 decades after the tran- 
sient; and (3) there is no dust and the star now has a very high photospheric temperature, ~ 70000 to 
100000 K, so that the faintness of the survivor is entirely due to bolometric corrections. In (1) and (2), the 
change in bolometric luminosity must be an order of magnitude or more. In case (3), such a radical increase 
in photospheric temperature would have be accompanied by significant mass loss which would be hard to 
reconcile with the requirement for no dust. Deeper ultraviolet observations than the available GALEX data 
would constrain this possibility, The remaining possibility is that the dust covering fraction is very small, 
less than 10%, with our line of sight coincidentally passing through one of the optically thick patches. The 
mid-IR emission is then reduced by the covering fraction. In this scenario we should still see the steady 
optical brightening created by an expanding shell. None of these possibilities seems terribly attractive or 
plausible. 
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If SN 1961V was a supernova, then it becomes one of the rare SNe with ob servations of its progenitor 
star, and the only one with a relatively detailed pre-explosion light curve (see ISmarttll2009l). We know it 



Goodrich et al 



was ve ry luminous, L* ~ 10 6 3 L Q , and likely hot in qui escence following the arguments of 
(119890 . Based on the emission line ratios reported by iGoodrich et al] dl989l) for the Western HII region 
(the Eastern HII region overlaps the region with SN 1961V and shows evidenc e for contamination by a 
SN remnant), we estimate an oxygen abundance of approximately 8.3 following iKewley & Ellisonl ([2008) 
or approximately ~ 1/3 Solar and similar to the metallicity of the LMC. This loc al measurement is a li ttle 
lower than estimates of ~ 1/2 Solar from the metallicity gradient measurements by lFerguson et all (11998b . If 
we select stars at the end points of the Padua (M arigo et al. isochrones with 10 6 A L Q < L* < 10 6 5 and 
20000 K < < 40000 K, they correspond to very massive stars, with Mzams > 80M Q for all metallicities 
from LMC to Solar. 

These pr operties are very simi l ar to those of the only othe r high mass progenitor to be identified, that 
of SN 2005gl (iGal-Yam et al.ll2007llGal-Yam & Leonardll20090 . For their measurements and parameters (a 
progenitor with V = 20.04 ±0.15 mag at 66 Mpc with E(B — V) ~ 0.07), this progenitor had a luminosity 
similar to that of SN 1961V, with U ^ 10 6 7 L for ?; = 20000 K. Like SN 1 961 V it was a Type Iln, and it had 
a comparable peak luminosity, near My ~ -17 mag. IGal-Yam et al.1 (120071) and lGal-Yam & Leonard! ([2009) 
propose that the spectral properties of SN 2005gl are best explained by heavy mass loss or mass ejections 
closely correlated with the SN. In fact, there is growing evidence that pre-supernova bursts of mass loss, 
while not common, are also no t rare. The most remar kable case is the eruption observed two years prior to 
the peculiar Type lb SN 2006jc (IPastorello et al. lboo7h . but the light curves of other SN show strong evidence 
for mass ejection episo d es shortly before t he SN (e.g., SN 2006gy, SN 2005ap, SN 2006tf, SN 2007va, see, 



e.g.. iSmith et all (|2008h . iKozlowski et all d2010h). The case of SN 2006jc may be particularly apt since the 
spectral evidence for excess helium ([Branch & Greensteinlll97l[) suggests that SN 1961V was close to being 
a Type lb rather than a Type Iln/pec. 

Suppose we interpret SN 1961V in this context. In this view, the pre-SN light curve of SN 1961V is 
mapped into the pre-SN history of mass loss, and thes e phases of mass loss are then mapped into the post- 
SN light curve. We modify the IGoodrich et all <|1989h scenario as follows. In quiescence, the progenitor 
is a hot star with a low density fast wind. Sometime before 1930 (likely closer to 1800), the star transi- 
tions from the compact, hot (T* ~ 40000 K) state with a low density, high velocity wind to (on average) 
a cooler (I 1 * ~ 7500 K) star with a high density, lower velocity wind. Then, around 1955 (not 1960) the 
star undergoes an LBV (or other) eruption to produce the pre-peak luminosit y plateau, accompanied by a 
furth er rise in the wind densit y and a higher wind velocity. The light curves in IB ranch & Greensteinl dl97ll) 
and lDoggett & Branchl d 19851) appear to allow the outburst phase to commence earlier than 1960 due to the 



poor quality of the magnitude limits from 1955 to 1960. Then in December 1961, the star undergoes core 
collapse and produces an SN, leading to the luminosity peak. The outgoing shock wave now interacts with 
the previous mass loss history, where the first, more luminous post-SN plateau is due to interactions with 
the LBV eruption ejecta, and the second, longer plateau is due to the wind emitted in the cool phase. The 
luminosity then drops dramatically when the shock wave reaches the low density wind of the pre-explosion 
hot star phase. 
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The observed velocities now represent the velocity of the expanding shock wave rather than the wind. 
Suppose we try to power the light curve using the luminosity available from shock heating the circumsteilar 
medium (CSM), 
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wher e M is the mass loss rate, v w is the wind velocity, v s is the shock velocity and e is the radiative efficiency 
(e.g. IChugai & Danziger 1994). Here we have scaled v s to the 4000 km/s suggested by the VLBI obser- 



Chu et al 



(12004I) . Using 2000 km/s simply drives the required mass loss rates upwards while 



vations of 

making it easier to have extended, post-SN luminosity plateaus. 

The difficult part about producing the first post-transient luminosity plateau is its duration. The length 
of the plateau t p should be roughly t p ~ t e v w /v s , where t e is the durati on of the eruption prio r to 1961. If 
t e ~ 1 year and t p ~ 0.5 years, as in the description of the light curve by lGoodrich et all ((1989), the required 
velocity ratio v w /v s ~ 1/2 seems unphysical. However, if however, the eruption commenced closer to 
1955 , so that t e ~ 5 years, but w as missed due to the shallowness of the observations (Branch & Greenstein 
197ll . lDoggett & Branchlll985r) . then we need only have v s /v s ~ 1/10. Suppose we adopt v s /v w = 5, then 
reproducing the plateau luminosity of order L ~ 2 x 1O 7 L (m pg ~ 17) requires a mass loss rate of 



M ~ 0.03 



/4000km/s\' 



v. 



■ 



(8) 



which is grossly consistent with an LBV eruption ( Humphreys & Davidson! 1994 ). This rises to lO _1 M0/year 
if we use v s = 2000 km/s. I n either case, enough m ass is involved that we may be also be underestimating 
the radiative efficiency (see lSmith & McCrayil2007r) . 



The shock then moves out through the lower density material from the earlier "S Doradus" phase - 
here we simply envision an extended period where the star is on average producing a relatively dense wind. 
The necessary shock luminosity is now 10 times lower, and we are free to make the ratio v s /v w much larger 
since we have no definitive time scale for the start of this phase beyond that it began before ~ 1930. Thus, 
the mass loss rate need only be 



M ~ 4 x 10" 



/4000km/s\ 



J VIOOkm/s 



0.1 



M /year, 



(9) 



which is relatively easy for an LBV produce even outside of eruptions. Here v s /v w = 40, so the enhanced 
mass loss phase would have started circa 1 800 in order to make t p ~ 4 years. The biggest problem with this 
schematic is that the material ejected prior to the SN cannot itself form significant amount s of dust or the 
progenitor would have been self-obscured, similar to SN 2008S (see, e.g.. |Prieto et al.ll2008l) . 



Not only was the progenitor of SN 1961V massive, Mzams <; 80M(D, but it must also h ave been rel- 
atively massive at death. While appearing to be rich in helium (IBranch & Greenstein! 1 1 97 II) . it was still 
a Type II SN, rather than a Type lb or Ic. We found no detai led pre-supernova models for this mass and 
metallicity range, but it is in the regime that lHeger et all (|2003) estimate would be weak Type Ib/c fall-back 
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SN leading to black hole formation. The solar metallicity models in lWoosley et al.1 (120021) have already lost 
much of their helium to mass loss. For their low metallicity models (10 -4 solar!), the star would need to be 
more massive at death than M ~ 30M© in order to retain any hydrogen. If it was a fall-back SN forming 
a BH, SN 1961V was not notably sub-luminous. Alternatively, some massive stars may still form NSs, 
as suggested by the existence of a magnetar in Westerlund 1. The progenitor of this NS seems to require a 
> 40M Q progenitor given the other ma ssive stars in the cluster (IMuno et al 1I2006L unless it can be explained 
by binary evolution and mass transfer ( Belczynski & Taam 20081) . Since it takes virtually no mass to power 
accretion onto a ~ 10M Q black hole at the Eddington limit compared to the ejected mass, we might expect 



the ne wly formed BH to accrete at the Eddington limit for an extended period of time. However. IPerna et al. 



(2008) found no X-ray emission from the site to a limit of Lx < 6 x 10 ergs/s (2-10kevfl in March 2000, 
corresponding to a limit of order 5% of Eddington. 

Finally, the existence of other Type Iln supernova requiring major mass loss events shortly before 
collapse changes our prior on the likelihood of such correlations for SN 1961V. Rather than being bizarre, 
it is simply the closest example, one which is so close that we could see the pre-SN activity. SN 1961V, 
SN 2005gl and their relatives are all cases where the correlated mass loss is large and dramatic. We usually 
assume that stars are evolving quasi-statically in their last phases, with no indicators of imminent death, yet 
this clearly does not hold for this class of objects. It is an interesting question whether this phenomenon 
is limited to a special class of SN, as proposed by Gal -Yam et all (120070 . or that we presently only notice 
the most dramatic examples of a more ubiquitous phenomenon. In either case, it appears that studies of SN 
progenitors should evolve from simple attempts to obtain a single snapshot of the star to monitoring their 
behavior over their final years. 
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With a similar limit of < 2 x 10 ergs/s for the softer 0.3-8 keV band (Soria & Perna 2010, private communication). 
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Table 1. Spitzer Observations of NGC 1058 



date 


MJD 


Pi/Program 


3.6/im 


4.5/im 


5.8/mi 


8.0/im 


24^m 


70/um 


2004-08-14 


53231.31 


Fazio/69 


150 


150 


150 


150 








2004-08-25 


53242.03 


Fazio/69 














10 


9 


2007-09-14 


54357.66 


Kotak/40619 


300 


300 


300 


300 








2007-09-29 


54372.22 


Kotak/40619 














30 






Note. — Exposure times are in seconds. The IRAC frame times were 30 s in both observations. 
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Table 2. IRAC Photometry 



OUUiLC 


A nprh i i*c 
r\pCI LUIC 


fii 




rs si 


rs m 






fmTvl 

\ mj yj 


I'm Tv "1 




\ mj y) 


SN1961V area 


I'.'A 


0.0139±0.0011 


0.0104 ±0.0014 


0.0557 ±0.0062 


0.1682 ± 0.0096 




3"6 


0.0226 ±0.0010 


0.0160 ±0.0019 


0.0690 ± 0.0074 


0.2070 ±0.0108 


SN1961V/#8 


DAOPHOT 


0.0078 ±0.0010 


0.0079 ±0.0010 


0.0288 ± 0.0098 


0.0614 ±0.0275 


SN1961V/other 


DAOPHOT 


0.0084 ±0.0008 


0.0081 ±0.0010 


0.0096 ± 0.0041 


0.0571 ±0.0064 


Star #3 


2'M 


0.0071 ±0.0013 


0.0031 ±0.0012 


< 0.019 


< 0.028 




3'.' 6 


0.0044 ±0.0021 


0.0003 ±0.0017 


0.0102 ±0.0084 


0.0349 ±0.0103 




DAOPHOT 


0.0103 ±0.0008 


0.0073 ±0.0010 






Star A 


2"4 


0.0190 ±0.0010 


0.0115 ±0.0010 


0.0346 ± 0.0070 


0.0238 ± 0.0078 




3'.' 6 


0.0204 ±0.0017 


0.0129 ±0.0016 


0.0566 ± 0.0075 


0.0167 ±0.0080 




DAOPHOT 


0.0256 ±0.0011 


0.0151 ±0.0013 






StarB 


2'.'4 


0.0147 ±0.0018 


0.0075 ±0.0014 


0.0209 ± 0.0092 


0.0529 ±0.0087 




3'.' 6 


0.0196 ±0.0029 


0.0119±0.0019 


0.0605 ±0.0081 


0.0869 ±0.0088 




DAOPHOT 


0.0149 ±0.0012 


0.0094 ±0.0011 






StarC 


2"4 


0.0064 ±0.0012 


0.0057 ±0.0011 


0.0367 ± 0.0051 


0.1031 ±0.0099 




3'.' 6 


0.0073 ±0.0019 


0.0041 ±0.0017 


0.0251 ±0.0063 


0.1273 ±0.0101 




DAOPHOT 


0.0056 ±0.0009 


0.0062 ±0.0012 






SN 2007gr 


2'M 


4.7079 ± 0.0207 


3.7626 ±0.0137 


2.7023 ±0.0431 


2.4228 ±0.1452 




3'.' 6 


5.1208 ±0.0313 


4.0373 ±0.0192 


3.0947 ±0.0225 


3.4095 ±0.1182 




DAOPHOT 


5.6602 ±0.1845 


4.1676±0.3363 






SN 1969L 


2'M 


0.0065 ±0.0013 


< 0.0060 


0.0093 ±0.0060 


0.0165 ±0.0101 




3'.' 6 


0.0099 ±0.0012 


< 0.0065 


0.0208 ± 0.0054 


0.0013 ± 0.0047 



Note. 



— Flux limits are 3a limits. 
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Table 3. MIPS Photometry 



Source 


[24] 


[70] 




(mJy) 


(mJy) 


SN1961V area 


0.226 ±0.0039 


<8.0 


Star #3 


0.089 ±0.0044 




Star A 


0.043 ±0.0082 




Star B 


0.094 ±0.0080 


< 190 


StarC 


0.027 ±0.0049 


< 178 


SN 2007gr 


3. 155 ±0.0398 




SN 1969L 


< 0.063 





Note. 



— Flux limits are 3a limits. 



